Dynamics of Overhauser Field under nuclear spin diffusion in a quantum dot 
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The coherence of electron spin can be significantly enhanced by locking the Overhauser field from 
nuclear spins using the nuclear spin preparation. We propose a theoretical model to calculate the 
long time dynamics of the Overhauser field under intrinsic nuclear spin diffusion in a quantum dot. 
We obtain a simplified diffusion equation that can be numerically solved and show quantitatively 
how the Knight shift and the electron-mediated nuclear spin fiip-flop affect the nuclear spin diffusion. 
The results explain several recent experimental observations, where the decay time of Overhauser 
field is measured under different configurations, including variation of the external magnetic field, 
the electron spin configuration in a double dot, and the initial nuclear spin polarization rate. 
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I. INTRODUCTION 

Electron spins in single quantum dots consist of one of 
the promising systems for realization of quantum com- 
putation yy. The spin state of a single electron in a 
quantum dot can be coherently controlled either opti- 
cally through fast laser pulses or electronically through 
tuning of gate voltages [U [5J [21 HI |S] . In experiments, 
the coherence time of the electron spin is limited by its 
hyperfine coupling to the nuclear spin environment in 
the host semiconductor material. The coupling causes 
spectral diffusion and gives a typical spin decoherence 
time T| ~ 15ns for the electron spin qubit [Sl[7]. This 
coherence time could be significantly prolonged with ap- 
plication of spin echo or other dynamic decoupling tech- 
niques [51. However, implementation of these techniques 
requires repeated applications of many laser pulses. Each 
pulse inevitably induces some noise by itself, which lim- 
its the practical performance of the dynamic decoupling 
techniques to suppress spin noise under real environ- 
ments. 

Another technique to increase the coherence time for 
the electron spin is through the dynamic nuclear-spin 
preparation (DNP) [H UHl El [HI [131 [Ml [15], which 
prepares the nuclear spin environment into certain con- 
figurations. The nuclear spins in this configuration col- 
lectively generate an effective magnetic field (the Over- 
hauser field) on the electron spin with small fluctuation, 
significantly reducing the spectral difi^usion of the elec- 
tron spin qubit. Although nuclear spins can be polarized 
through many methods such as optical pumping [5] , sub- 
stantial reduction in the fluctuation of the Overhauser 
field requires almost complete polarization of the nuclear 
spins [TD], which is hard to achieve experimentally. The 
recent experiments, however, demonstrate a surprising 
feedback mechanism which can lock the Overhauser field 
to certain values without significant polarization of the 
nuclear spin environment [T^ [T31 [TH [T5]. The Over- 
hauser field generated from this locking mechanism has 
small fiuctuation, which effectively increase the coherence 
time T2 of the electron spin qubit by up to two orders of 



magnitudes. The locking of the Overhauser field through 
the DNP process has been seen in a number of experi- 
mental systems, under both optical and electronic gate 
control [laiiallllllll]. 

To enjoy a longer spin coherence time, it is desir- 
able to do gate experiments on the electron spins under 
the nuclear spin environment prepared with the DNP, 
which produces a fixed Overhauser field. An important 
question is then how long this fixed Overhauser field 
can survive after the DNP process and what facts de- 
termine/influence the relaxation time of the Overhauser 
field. Recent experiments have observed that the pro- 
duced Overhauser field can endure from a few seconds to 
a few minutes and the relaxation time of the Overhauser 
field depends on a variety of experimental parameters, 
such as the applied magnetic field {T6\ [13 [ISj , the elec- 
tron spin configuration in double quantum dots |17j , and 
the DNP pump time . 

In this paper, we develop a quantitative theory to cal- 
culate the relaxation time of the Overhauser field un- 
der the environment of quantum dots, and provide an 
explanation to the experimental observations mentioned 
above in different setups and with control of different 
experimental parameters. The relaxation of the Over- 
hauser field is caused by nuclear spin diffusion, which 
has been well studied in bulk material [201 HI], where 
diffusion is caused simply by the nuclear dipole-dipole 
interaction. In a quantum dot, however, the presence 
of the electron spin generates several new effects. First, 
the electron spin can mediate the diffusion of the nuclear 
spins through a virtual hyperfine process. Second, the ef- 
fective magnetic field from the electron spin produces an 
inhomogeneous Knight shift on the nuclear spins through 
the hyperfine coupling which suppresses the nuclear spin 
diffusion. The change of the nuclear spin diffusion coef- 
ficient by the Knight shift has been taken into account 
in Ref. [25], but without consideration of the electron 
mediated nuclear spin diffusion. A recent work considers 
relaxation of the Overhauser field due to the electron me- 
diated nuclear spin diffusion, but without consideration 
of the direct nuclear dipole-dipole interaction [53^ . In this 



2 



case, the Overhauser field can only decay by less than 1%, 
which suggests that it is necessary to include the dipole 
interaction in the long time dynamics. A quantitative 
theory is still missing to our knowledge that includes a 
complete description of all the competing effects men- 
tioned above. In this work, we take into account all these 
diffusion mechanisms, and the resulting theory provides 
an explanation of the recent experimental observations 
inRefs. jlS HH HH US] . 

The paper is arranged as follows: in Sec. II we give 
a formalism to describe relaxation of the Overhauser 
field that includes contributions from the nuclear dipole- 
dipole interaction, the electron mediated nuclear spin dif- 
fusion, and the Knight shift. The effective nuclear spin 
diffusion equation is solved numerically to determine the 
relaxation time of the Overhauser field. In Sec. Ill, we 
compare the theoretical calculations with the recent ex- 
perimental observations under control of different experi- 
mental parameters, and show that they are in qualitative 
or semi-quantitative agreement. We summarize our re- 
sults in Sec. IV with brief discussions. 



II. DECAY OF THE OVERHAUSER FIELD 
THROUGH NUCLEAR SPIN DIFFUSION 

We assume that an external magnetic field Bq much 
larger than the mean value and variance of the local 
Overhauser field generated by nuclei is applied along the 
z-direction (perpendicular to the quantum dot layer) . In 
this case, we can drop the nonsecular terms in the inter- 
action Hamiltonian |24j . For simplification, we consider 
only one species of nuclei around the quantum dot elec- 
tron. The total Hamiltonian for the electron and nuclear 
spin system, including both the Fermi contact hyperfine 
interaction and nuclear dipole-dipole interaction, can be 
written as: 



H = He + Hn + Hen + Hnn, 
Hn = —gnf^NBp Ij, 



(1) 

(2) 
(3) 



(5) 



B^j = g(g„Miv)'i?,y'(l-3cos2%), (6) 

where Ai denotes the hyperfine coupling between the 
electron and nuclear spin at site i with spatial coordi- 
nates {xi,yi,Zi). Rij is the distance between two nuclei 
at site Oij is the angle between the line connecting 
sites i,j and the z direction. 

We note that for Bq ranging from a few mT to a 
few T, the electron Zeeman splitting is on the order of 



10~^ — lO^GHz, while the average hyperfine coupling in 
most quantum dot systems is on the order of MHz. Thus 
we can adiabatically eliminate the spin-flip terms in the 
hyperfine interaction Hamiltonian and correspondingly 
modify the other terms in the Hamiltonian as [23j: 
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where we have introduced an electron-mediated nuclear 
flip-flop term in i/„„. Since we are interested in the long 
time dynamics of nuclear spins, we can completely elimi- 
nate the electron from the Hamiltonian by replacing the 
constant operator with its expectation value. We find 
that using — 1/2 or — —1/2 will yield almost the 
same result in the following calculations. Therefore we 
can set — 1/2 for simplicity and arrive at the following 
effective Hamiltonian (neglecting constant terms): 

Ho + Hi, (10) 

i-gnmBo + A,/2)It + 2S,,/;//,(ll) 



H 
Ho 



E 



H, = -5](S,, 



AiAj 



(12) 



Here, the term proportional to Ai in Hq is the Knight 
shift term. For this Knight shift, we have neglected the 
small term proportional to Af in Eq. ([8| as it is domi- 
nated by the Ai term. 

The expectation value for z component of the nuclear 
spin at site k will evolve according to the Schrodinger 
equation: 



dm 



-Tr{p(i)[iJi,/f]}, 



(13) 



dt 

where p{t) is the nuclear spin density matrix at time t, 
which can be calculated by switching to the interaction 
picture: 

m = p{o) + y\m,H,{t')]dt', (14) 

with Hi{t) — exp{iHot/h)Hiexp{—iHot/h). Further cal- 
culation yields pUj : 

d{I, 



^ = -Tr{pmHiit),II]} 

'''Tr{p{t^t')[H,,[H,{t'),m]}dt' (15) 





We assume the nuclear spin (with spin-/) density ma- 
trix in a product state of the following form: 



p{t) = <^Pk{t) Pk{t) = 



1 



m)) 
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(16) 
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Such an approximation is valid for the initial state of nu- 
clear spin when correlations and transverse coherence are 
negligible in its equilibrium configuration. Off-diagonal 
terms may appear in p(t) for t > 0, but it can be shown 
that they only have a minor contribution to the evolution 
of {I^{t)) compared to the diagonal part expressed in Eq. 
(fTel (see Ref. HO])- 



expansion of {I^{t)) for site i around site k: 



2 dx-dxP (^fe-^JK-^j- 



(22) 
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By using the explicit form of the Hamiltonian [Eq. ( 10 
and density matrix [Eq. ([16])], we can reduce Eq 



where Einstein's summation convention is implied for 
spatial index ck,/3. Substituting this into Eq. (17) and 



(151 to: 



noting that the summation of the first order derivative 
term over all sites vanishes due to the lattice symmetry 
1251. we have: 



dt 



(17) d{Il) 



dt 



(23) 



Tr{ [Hi (t) , /|] [Hi (t ~ t'), If] }dt' The notation above means summarization over the 

sites near k. Define coefficient D"^ = ~ 



where Wki has a clear physical meaning as the ffip-ffop 
rate between nuclear spins at site i and k. 

For a 2D InAs/GaAs quantum dot, we take As nuclei 
(/ — 3/2) as an example for further calculation. The 
parameter Wki can be analytically calculated when ap- 
proximating the integration upper limit in the above ex- 
pression for Wki from t to infinity |22] . 
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The hyperfine coupling rate Ai is proportional to square 
of the electron wave function in a quantum dot. And in 
the following calculation we assume the dot potential is 
like a square well in the z-direction and the electron wave 
function takes a Gaussian shape in the x, y-plane. The 
hyperfine coupling rate Ai can then be written as Ai = 
Aocos^(7rzi/zo) exp[-(a;f -h yD/Zo], where ( ) are 

spatial coordinate of the site i. Iq and zq are, respectively, 
the Fock-Darwin radius and thickness of the quantum 
dot, and Aq is the hyperfine coupling for the nuclear spin 
at the origin where the electron locates. 

We assume the nuclei spins follow a diffusion process, 
where the ffip-ffop rate for two distant sites is negligible. 
This approximation is justiffed by the fact that the coeffi- 
cient Cik generally decays to zero fast as Rik increases. If 
we treat {I kit)) as a continuous function of spatial vari- 
able a;" (a = x,y,z). We can then carry out a Taylor 



)/2 and similarly note that for a f3 the 



summation over all sites vanishes, we have: 



dt 



dx'^ 
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dz- 



){im) (24) 



Eq. (24 1 is a 3D anisotropic diffusion equation with 



spatially varying diffusion coefficients (as Aik, Bik, Wik 
all depend on the spatial coordinates), which is not easy 
to solve. To further simplify it, we note that to obtain 
the major feature for the full time dynamics of the Over- 
hauser field {hz{t)) = J2k ^k{Ik{t)) ^ it is reasonable to 
first ignore the diffusion in the z direction as the quan- 
tum dot layer is usually a few nm thick and chemical 
or structural mismatch in adjacent layers may strongly 
suppress diffusion in the z direction [18 . In addition, 
from symmetry in the 2D x, y-plane, we expect to have 
jjxx ^ jjyy thus define an average 2D diffusion co- 
eflicient D{x,y) = J2zKk Wik[{xk - x,)^ + {yk - Vif]/^- 
Now we have a simplified 2D diffusion equation: 



dt 



D{x,y){ 



d^ 



d- 



dx^ dy 



(25) 



The above partial differential equation can be effective 
solved using the finite element method by coarse graining 
a large number of nuclear spin sites to a small number 
of mesh nodes. But before solving Eq. ( 25 1 , we would 



like to have some discussion about the diffusion coeffi- 
cient D{x,y). For x,y ^ Iq, the role of electron can be 
neglected and a numerical calculation of the above diffu- 
sion coefficient gives a uniform background value D 7 
nm? / s, which is consistent with the previous theoretical 
and experimental study of the diffusion coefficient in the 
bulk material [iniHIlHl]- In our calculation, we set the 
quantum dot parameters as Iq = 30 nm, zq = 10 nm, 
Ao = 1 ^leV « 1.5 MHz, Efe^fc ~ 80 ^leV, the lat- 
tice constant oq — 0.563 nm, and the number of nuclei 
N ^ 9x 10^, according to the typical experimental values 
[16lll7l[l8]. 

Within the range of Fock-Darwin radius Iq, however, 
the presence of the quantum dot electron will change the 
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diffusion coefficient through two competing mechanisms: 
on the one hand, the presence of the confined electron 
generates an inhomogeneous Knight shift [3], which lifts 
the degeneracy of the nuclear Zeeman energy for differ- 
ent nuclei and prevents the spin flip-flop; on the other 
hand, electron mediated nuclear spin flip-flop enhances 
the nuclear spin diffusion from the center to the edge of 
the dot. 

Our numerical simulation shows that whether one 
mechanism dominates the other is largely determined by 
the external magnetic field Bq. Fig. [T][2]show the diffu- 
sion coefficient D{x,y) under Bq = 0.2T and Bq — 2T. 
We can see that the electron mediated flip-flop greatly 
enhances the nuclear spin diffusion near the center of the 
dot under a small magnetic field (with a sharp peak in 
D{x,y)), while for a large magnetic field, this effect is 
negligible compared to the Knight shift which suppresses 
the nuclear spin diffusion (with a wide dip in D{x,y)). 
The difference can be easily explained from the effective 
Hamiltonian [Eq. ( 10]|T2| ]: the electron mediated flip- 
flop term is inversely proportional to Bq while the Knight 
shift term is independent of Bq. We note that the rea- 
son why we have a narrower peak than the dip is due to 
the fact that the Knight shift term is proportional to the 
hyperfine coupling (with a Gaussian distribution) while 
the Electron mediated flip-flop term is proportional to 
the product of two nuclei's hyperfine coupling rates. We 
also note that the 2D diffusion coefficient D{x,y) in Fig. 
[T]|2]does not have azimuthal symmetry since in our cal- 
culation we assume the nuclear spins in a square lattice 
which has no azimuthal symmetry. 
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FIG. 1: Diffusion Coefficient D(x,y) under Bo = 0.2T. The 
narrow high peak at the center of the dot is due to electron- 
mediated nuclear spin flip-flop, and tlie wide dip is due to 
inhomogeneous Knight shift. 
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III. COMPARISON WITH EXPERIMENTS 

To compare with experiments, we numerically solve 
the diffusion equation [Eq. (25 1] under certain initial 
and boundary conditions. For the initial condition, the 
nuclear spins are partially polarized through the DNP 
process from the hyperfine interaction with the elec- 
tron spin [T^ [T3] [T3J [TS]. It is reasonable to expect 
that right after the DNP process, the polarization dis- 
tribution (/^) is proportional to the hyperfine interac- 
tion rate. So in the following calculation, we assume 
(J^) cx Ak (X exp[-(x^ -I- y])/ll] at i = for solving 



the diffusion equation [Eq. (25l]. We can take the natu- 



ral boundary condition where (/|) approaches zero when 
the radius goes to infinity. However, in numerical calcu- 
lation, we have to take a finite radius. To make the spin 
diffusion possible, this finite radius has to be significantly 
larger than the radius of the size of the initial electron 
wave packet. In the calculation, we typically take the 
integration size about 300 nm (10 times of Iq) so that 
the number of total nuclear spins in this area is about 
100 times the number of initially partially polarized nu- 
clear spins covered by the electron wave packet. With 
this choice, we expect the cutoff error should be small 
(at a percent level) . 
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FIG. 2: Diffusion Coeflicient D(x,y) under Bo = 2T. The in- 
homogeneous Knight shift dominates in this case, so diffusion 
is generally suppressed within the Fock-Darwin radius. 



First, to compare with the experiments in Ref. |17j . we 
calculate the relaxation of the Overhauser field h^{t) = 
^k{Ik{t)) under different electron states and different 
values of the external magnetic field Bq, and the result 
is shown in Fig. [3) We note that for the double quan- 
tum dot system in Ref. [T7] , if the electron stays in the 
(2,0) singlet state, the electron spin has Sz = 0, with 
basically no influence on the nuclear spin diffusion. In 
this case, the nuclear spin diffusion is governed by the 
intrinsic nuclear dipole-dipolc interaction. However, for 
the electron in the (1,1) state, with the magnetic field 
in the range of tens of mT as in this experiment, the 
electron mediated spin diffusion dominates the Knight 
shift and it accelerates the nuclear spin relaxation. One 
can see two effects from Fig. [3] with the magnetic field 
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FIG. 3: Decay of the Overhauser field under various small 
magnetic fields. The solid line refers to the case with electron 
staying in the (2,0) singlet state where electron plays no role 
in nuclear spin diffusion. 



FIG. 5: Overhauser Decay of the Overhauser field under vari- 
ous initial distributions of the nuclear polarization. The solid 
line refers to the case with a long DNP pump time that gives 
a Gaussian distribution with the size characterized by the 
Fock-Darwin radius lo- Other lines correspond to narrower 
polarization distribution characterized by a Gaussian with its 
size ro < lo- 
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FIG. 4: Decay of the Overhauser field under a strong magnetic 
field. The solid line shows the dynamics of the Overhauser 
field under B — 2T and the dotted line is a fit by using a 
constant diffusion coefficient. 



Bo = 10 mT or 20 mT: (1) electrons in the (1,1) state will 
speed up the decay of the Overhauser field compared to 
electrons in the (2,0) state; (2) a smaller magnetic field 
gives a faster decay of the Overhauser field. Both of these 
effects agree with the experimental observations in Ref. 
[17j , and the decay time scale is also consistent with what 
one measures from the experiments in term of the order 
of magnitude. 

With a much larger magnetic field (say, Bq = 2T, as in 
experiments in Ref. |18j). the electron mediated nuclear 
spin diffusion is suppressed, and the Knight shift plays a 
more important role. The Knight shift suppresses the nu- 
clear spin diffusion, and can prolong the relaxation time 



of the Overhauser field to make it significantly larger than 
the relaxation time in the bulk material. Fig. [4]shows de- 
cay of the Overhauser field in this case, and we can fit the 
curve with an effective constant diffusion coefficient with 
its value about -De// ~ 0.7nm^ /s. Compared with the 
diffusion coefficient in the bulk material (D Inm? /s), 
the effective diffusion here is suppressed by a factor of 10 
under a strong magnetic field Bq. Experiments done in 
Ref. [TB] yield an effective diffusion coefficient 50 times 
smaller than the value in the bulk system. While the 
suppression there could have contribution from other fac- 
tors, such as inhomogeneity of lattice constant there, we 
believe that the Knight-shift induced suppression plays 
a large role in this experiment under a strong magnetic 
field. 

The experiment in Ref. [16] studies relaxation of the 
Overhauser field under different pumping time for the 
DNP process. With a shorter DNP pumping time, the 
nuclear spin polarization may have a narrower distribu- 
tion in space [16]. Although we do not know the exact 
distribution of the nuclear spin polarization from a short 
DNP pump, to model this effect qualitatively, we simply 
assume that this distribution (/|) is still a Gaussian but 
with its radius <Iq. Taking this (/|) as the initial con- 
dition, we can calculate relaxation of the corresponding 
Overhauser field from the diffusion equation [Eq. (25l], 
and the result is shown in Fig. [5] The result indicates 
that a narrower distribution of initial nuclear spin polar- 
ization leads to a faster decay of the Overhauser field, 
which is consistent with the experimental result in Ref. 
[16] . This effect can be explained by noting that the dif- 
fusion is much stronger near the center of the dot due to 
the electron mediated diffusion peak (see Fig. [ij, which 
results in a shorter relaxation time of the Overhauser 
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field if the initial polarization is more concentrated near 
the dot center. 



IV. SUMMARY AND DISCUSSION 

In summary, we have established an effective method 
for calculating the long time dynamics of the Overhauser 
field under nuclear spin diffusion and shown that the con- 
fined electron in a quantum dot can both enhance decay 
of the Overhauser field by mediating nuclear spin flip-flop 
and suppress nuclear spin diffusion via inhomogeneous 
Knight shift. Which effect dominates depends critically 
on the magnitude of the external magnetic field. We sim- 
ulate the relaxation of Overhauser field under different 
electron spin configuration, external magnetic field, and 
initial nuclear polarization distribution, and the results 
agree qualitatively with a series of recent experimental 
observations. 

For the purpose of maintaining the Overhauser field 
generated by the DNP process as long as possible, ap- 
plying a large magnetic field turns out to be the most 
effective way since a large magnetic field suppresses the 
electron mediated nuclear spin diffusion. The intrinsic 
nuclear spin diffusion from the nuclear dipole-dipole in- 
teraction is also suppressed by the inhomogeneous Knight 
shift. These two kinds of suppression, combined together. 



leads to a long relaxation time of the Overhauser field. 

For calculations in this paper, we focus on the relax- 
ation dynamics of the expectation value of the Over- 
hauser field, since this is the quantity that has been 
measured in several recent experiments. Similar meth- 
ods could also apply to calculation of relaxation of the 
variance of the Overhauser field, and in term of the time 
scale, the relaxation time for variance is basically the 
same as the relaxation time for the expectation value. 

During the DNP process, nuclear spin diffusion also 
takes place. The final distribution of the nuclear spin po- 
larization and the limit on its fluctuation may depend on 
a balance between the DNP pumping cycle and the con- 
tinuous nuclear spin diffusion process |T5J [121 IHl ESI HZ] ■ 
To understand this balance, we need a more detailed un- 
derstanding of the dynamics for the DNP process and its 
dependence on various experimental control parameters. 
This is an interesting topic for further investigation. 
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